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*> : 1 Introduction 

Loop corrections modify the Coulomb potential: electron loop insertion into the photon 
1 propagator leads to the Uehling-Serber correction to the electric potential of point-like nuclei 
Q\- pp. Though it leads to an important phenomenon contributing to the Lamb shift of the 
energies of atomic electrons numerically the shift being of the order of a 5 m e is small (we 
are using Gauss system of units, where a = e 2 = 1/137 and in all formulas h = c = 1 is 
implied). 

Analogous correction in the case of external magnetic field qualitatively change the be- 
havior of atomic energies: in particular the energy of the ground level remains finite in the 
limit B — > oo; also spontaneous production of positrons becomes possible only for nuclei with 
Z > 52. Without taking radiative corrections into account in the limit of infinite magnetic 
field energy of ground atomic level tends to minus infinity and point-like nucleus with any Z 
becomes critical at large enough B. At magnetic fields B > (Ze 2 ) 2 m 2 /e the characteristic size 
of the electron wave function in the transverse to the magnetic field direction an = 1/ V eB 
(the so-called Landau radius) becomes smaller than Bohr radius a# = l/(Ze 2 m e ) making 
the Coulomb problem essentially one-dimensional. Singularity of the Coulomb potential in 
d — 1 is stronger than in d = 3. In d = 1 the energy of the ground level is unbounded from 
below: the "fall to the center" phenomenon occurs. In the case of external magnetic field the 
singularity is cured by the finite value of a#: at \z\ < a H the Coulomb problem remains three 
dimensional. This is the reason why ground level goes down when B grows. At superstrong 
magnetic fields B > "iTxm 2 /e z radiative corrections screen the Coulomb potential at short 
distances \z\ < l/m e and the freezing of a ground state energy occurs: it remains finite at 
B — > oo ([2], 13], PQ). For Z > 52 the value of freezing energy is below — m e , so the ground 
level enters lower continuum when B increases and spontaneous production of e + e~ pair 
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from vacuum becomes energetically possible and thus takes place. In this process electron 
occupies ground level while positron is emitted to infinity. For Z < 52 freezing energy is 
above — m e and spontaneous positron production does not occur [5]. 

There exists the direct correspondence between radiative corrections to the Coulomb 
potential in d = 3 case in strong magnetic field B > B = m 2 e /e and radiative corrections 
to the Coulomb potential in d = 1 QED. That is why we start our presentation (in Section 
2) from the analysis of the Coulomb potential in D = 2 QED of massive fermions. When 
these fermions are light, g 2 > m 2 , the exponential screening of the Coulomb potential at 
short distances occurs. In the limit m — > (massless D = 2 QED, the so-called Schwinger 
model) this exponential screening occurs at all distances because photon gets mass, m 7 = 2g 
[5J. In Section 3 we analyze radiative corrections to the Coulomb potential in D = 4 QED 
in external magnetic field. The role of the coupling constant g 2 here plays the product 
e 3 B, and for e 3 B > m 2 the screening of the Coulomb potential occurs as well. In Section 
4 the structure of atomic levels on which the lowest Landau level (LLL) in the presence of 
atomic nucleus splits is determined. In Section 5 the Dirac equation for hydrogenlike ion at 
superstrong magnetic field will be derived and effect of screening will be studied for Z — 1. 
In Section 6 the influence of the screening of the Coulomb potential on the values of critical 
nuclei charges is discussed. In Section 7 the obtained results are summarized. 

Let us finish the Introduction discussing the numerical values of magnetic fields we are 
dealing with in these lectures. The magnetic field at which the Bohr radius of a hydrogen 
atom becomes equal to Landau radius is B a = e 3 m 2 rs 2 • 10 9 gauss, which is much larger 
than a magnetic field ever made artificially on Earth: B^ w 3 ■ 10 7 gauss. An interest to 
the atomic spectrum in the magnetic fields B > B a was triggered by the experiments with 
semiconductors, where electron-hole bound system called exciton is formed. Both effective 
charge and mass of electrons in semiconductors are much lower than in vacuum making B a 
in kilogauss scale reachable. 

The so-called Schwinger magnetic field Bq = m 2 /e pa 4.4 • 10 13 gauss and magnetic field 
at which the screening of the Coulomb potential occurs B pa 37rm 2 /e 3 pa 6 • 10 16 gauss 
should be compared with the magnetic fields at pulsars ~ 10 13 gauss and magnetars ~ 10 15 
gauss. Although the application of the results obtained in the condensed matter physics (say 
graphen, where the mass of charge carrier can be arbitrary low while the value of charge 
approach one) can not be excluded, our main interest in the problem considered is purely 
theoretical. 



2 The Coulomb potential in D = 2 QED of massive 
fermions 

Summing up diagrams shown in Fig. 1 we get the following formula for the potential of 
point-like charge: 

*(*) = Mk) = - fc2 ^f (fc 2) > n - - - ^Mk 2 ) , (i) 

where U(k 2 ) is the one-loop expression for the photon polarization operator. It can be 
obtained from the textbook expression for the polarization operator calculated with the help 
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of dimensional regularization [7j in the limit D — > 2: 

1 



m 2 ) = v 



t(i+t) 



ln(vTTI + Vt) - 1 



-4g 2 P(t) 



(2) 



where t = —k 2 /4m 2 , [g] =mass. 



^\/\/V 



+ 



+ 



Fig. 1. Modification of the Coulomb potential due to the dressing of the photon propagator. 

Let us note that U(k 2 ) is finite though corresponding integral ~ / d 2 p/p 2 is divergent in 
ultraviolet. The point is that the trace of gamma matrices which multiplies divergent integral 
is zero. In dimensional regularization the trace is proportional to D — 2 while ultraviolet 
divergency of integral over virtual momentum produces the factor 1/(D — 2) and the product 
of these two factors is finite. 

In order to obtain an expression for the Coulomb potential in the coordinate representa- 
tion we take k = (0, fcy) and make the Fourier transformation: 



ing 



e ik w z dk\\/2TT 



kj + 4g 2 P(k 2 /4m 2 ) 

The potential energy for the charges +g and — g is 

V(z) = -g$(z) . 



(3) 



(4) 



The integral in ([3]) cannot be expressed through elementary functions. However it appears 
possible to find an interpolating formula for P(t) which has good accuracy and is simple 
enough for the Fourier transformation to be performed analytically. 

The asymptotics of P(t) are: 



P(t) 



.3 ' 



t < 1 
t > 1 



Let us take as an interpolating formula for P(t) the following expression: 

2t 



P(t) 



(5) 



(6) 



3 + 2t 

The accuracy of this approximation is not worse than 10% for the whole interval oft variation, 
< t < oo. Substituting an interpolating formula in (j3]) we get: 



Airg 



e lk w z dk\\/2ii 



k\ + Ag 2 (kl /2m 2 )/(3 + fcjf /2m 2 ) 



4-rrg 



2g 2 /3m 2 



+ 



2g 2 /3m 2 



k\ k\ + 6m 2 + 4g 2 



Ak\\ z 



dk\\ 

~2^ 



(7) 



4-rrg 



1 + 2# 2 /3m 2 



1 

~2^' 



g 2 /3m 2 



^6m' 2 + Ag' 2 



=exp( — a/ 6m 2 + 4g 2 |,2|) 
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In the case of heavy fermions (m ^> g) 
the corrections are suppressed as g 2 /m 2 . 
In the case of light fermions (to -C g)- 



the potential is given by the tree level expression; 



*(*) 



z < - In 



z > i In ( 2- 

g 



ne ~2 9 \z\ 

m^g y -2ng (f^) \z\ 

m = corresponds to Schwinger model; photon gets mass. 

Light fermions make transition from m > g to m = continuous, 
fermions the Coulomb potential in D = 2 QED is screened at distances 

In Fig. 2 the potential energy for g = 0.5, m = 0.1 is shown. 
V(0) = 0. 



(8) 



In the case of light 

M £ V(2<7). 

It is normalized to 



V,V 




z 



Fig.2. Potential energy of the charges +g and —g in D = 2. The solid curve corresponds 

to P; the dashed curve corresponds to P. 

3 The Coulomb potential in D = 4 QED in superstrong 
magnetic field 

In order to find the potential of a point-like charge we need an expression for photon polar- 
ization operator in the external magnetic field B. Long ago an expression for the electron 
propagator in constant and homogeneous external magnetic field was found by Schwinger 
[S] as a parametric integral. For B ^> B = m 2 e /e integration can be easily performed and 
compact expression for G(k) follows. Using it one obtains an analytic expression for a photon 
polarization operator (see for example [9]). 

To understand the reason for great simplification of the expression for the electron prop- 
agator in the limit B B one should start from the spectral representation of the propaga- 
tor. The solutions of Dirac equation in the homogeneous constant in time B are known, so 
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one can write the spectral representation of the electron Green function. The denominators 
contain k 2 — m 2 — 2neB, and for B » m 2 /e and k 2 << eB in sum over levels the lowest 
Landau level (LLL, n = 0) dominates. In the coordinate representation the transverse part 
of LLL wave function is: \1/ ~ exp((—x 2 — y 2 )eB) which in the momentum representation 
gives \1/ ~ exp((—k 2 — k 2 )/eB) (we suppose that B is directed along the z axis). 

Substituting the electron Green functions we get the expression for the polarization op- 
erator in superstrong B. 

For B » Bq, k 2 << eB the following expression is valid [TO] : 



IU ~ e 2 eBl^eM # + Q « 



eB iy eB ' 

(q + k) 2 x + (q + k) 2 1 n . , fQ . 

* exp( -)dq dq z ^^- (1 - 27172)7^ * (9) 

eB q 0>z - m 

* - r (1 _ llll2) = e ' B * exp (-^i) * n 2^n = fc *) • 

Qo,z + k ,z -m leB 
With the help of it the following result was obtained in [4]: 

*W = 3 ^ , 2 x , (10) 

k 2 + k\ + ^exp P f ^ 



$(z) 



47re 

e 



e ifc ii z rfA;||d 2 A; ± /(27r) J 



fc J + k l + ^r^M-k 2 ± /{2eB)){k 2 /2m 2 e )/{3 + £f /2m| 

6^? e \z\ + e -yJ{2/-K)e 3 B+& m l\z\ 



1-e 



:n) 



For the magnetic fields B <C 37im 2 /e 3 the potential is Coulomb up to small power 
suppressed terms: 



e 3 B <C ml 



e 
\z\ 



1 + 



e 6 B 



(12) 



in full accordance with the D = 2 case with the substitution e 3 B — > g 2 . 

In the opposite case of the superstrong magnetic fields B 3> 3imi 2 /e 3 we get: 



( * e <rVW*&BM) 1 ]n(Jf%) >\z\>-^ 

1 ln(\/^) 



e 

1*1 



(1 _ e (-A/6mfkD) ? 



> 2 > 



(13) 



z > 



V^) = -e*(s) . (14) 

The Coulomb potential is screened at short distances l/m e > \z\ > 1/v ' e 3 B = a#/e. 
In Fig. 3 the plot of a modified by the superstrong magnetic field Coulomb potential as 
well as its short- and long-distance asymptotics are presented. 
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Fig. 3. The modified Coulomb potential at B = 10 17 G (blue, dark solid) and its long 
distance (green, pale solid) and short distance (red, dashed) asymptotics. 

Let us find the 3-dimensional shape of the screened Coulomb potential. The behavior 
of the potential in the transverse plane (z = 0) can be found analytically in the limit 
B 3> 3-Kml/e 3 from the general expression: 



$(z,p) = Aire J 



e i ^P +ik \\ z dk ]l d 2 k ± /{27r 



neglecting exponent in the denominator, which is valid for p > a#: 



(15) 



$(0,p) 



37rm| e 
e z B p 



^(2/n)e 3 B V 



< P 



(16) 



and the Coulomb potential is screened at large p in complete analogy with the D = 2 



case. 



For \z\ ^> l/m e in the integral (|15|) the values \k»\ <C m e dominate and we get: 

e 



:>l/m e 



z 2 + (i + €^) P 



Storm; ' 



(17) 
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Fig. 4. The equipotential lines at B/B = 10 4 . The dashed line corresponds to yjz 2 + p 2 



1 :ln 



y/{2/n)e3B V 3nmi " 

In Fig. 4 the equipotential lines are shown. The behavior of the screened Coulomb 
potential in the transverse plane was found numerically in [2]. 
Finally for 37rm 2 /e 3 3> B ^> m 2 /e expanding ( II 7p we get: 

$(r) = -(l-— (£AB )sin 2 fll , 6 = f A B, r = J p 2 + z 2 > l/m e , (18) 
r I 67T J v 

which coincides with the result obtained in [10] where the expression for the photon polar- 
ization operator at B > B was obtained as well. 

The expression for the screened Coulomb potential was obtained from the one-loop con- 
tribution to the photon polarization operator in the external magnetic field S > Bo. In 
momentum space it looks like: 

Aire 

$(k h k = k ± = 0) = . (19) 



'II vr fc^+6m,2 

,2 



If higher loops contain the terms ~ e 3 B (e 3 B / fc 2 ) n_1 they will drastically change the 
shape of the potential in the coordinate space. 

To calculate the radiative corrections one should use the electron propagator in an exter- 
nal homogeneous magnetic field B. The spectral representation of the electron propagator 
is a sum over Landau levels and for B 3> -Bo the contribution of the lowest level dominates 

mm 

G W = e -i/.B (1 _ i7l72) |i±^ _ (20) 

^0,3 " l e 
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where the projector (1—27172) selects the virtual electron state with its spin oriented opposite 
to the direction of the magnetic field B = (0, 0, B). The contributions of the excited Landau 
levels to G yield a term in the denominator proportional to eB and they produce a correction 
of order e 2 = a in the denominator of (|19p . 

Two kinds of terms contribute to the polarization operator at the two- loop level. First, 
there are terms in the electron propagators which represent the contributions of higher 
Landau levels. Just like in the one-loop case they produce corrections suppressed as e 2 
in the denominator of (TL9j) . i.e. terms of the order e 5 B which can be safely neglected in 
comparison with the leading ~ e 3 B term. Second, there is the contribution from the leading 
term in the electron propagator, given by (1201) . Let us consider the simplest diagram: the 
photon dressing of the electron propagator. Neglecting the electron mass we get: 

7^(1 - ^7172)^0,37^ = -2fe,3 - ^0,37271] = -2&o,3(l + 27172) , (21) 

which gives zero when multiplying external propagator of electron, since (1 + 27172) (1 — 
27172) — 0. This result is a manifestation of the following well-known fact: in massless QED 
in D — 2 (Schwinger model) all loop diagrams are zero except the one-loop term in the 
photon polarization operator (see for example p2]). That is why the contributions of the 
second kind are of the order of a(e 3 B) + m 2 ) 2 ) and they are not important. 

The generalization of the above arguments to higher loops is straightforward. Let us 
note that absence of higher loop corrections to polarization operator in Schwinger model 
is related to the absence of renormalization of axial anomaly by higher loops. In D = 2 
anomaly is given by correlator of two currents and axial current is proportional to vector 
current (see for example [T3"]). 

4 Atomic levels in superstrong B 

We are interested in the spectrum of a hydrogen-like ion in a very strong magnetic field B. 
We will write all formulas for hydrogen since their generalization for Z > 1 is straightforward. 

In the absence of magnetic field the spatial size of the wave function of the ground 
state atomic electron is characterized by the Bohr radius a# = l/(2n e e 2 ), its energy equals 
E = — m e e 4 /2 = — Ry, where Ry is the Rydberg constant. The transverse (with respect to 
B) size of the ground state of the electron wave function in an external magnetic field B is 
characterized by the Landau radius an — 1/ y/eB. The Larmour frequency of the electron 
precession is ujl = eB/m e . For a magnetic field B a = e 3 ml = 2.35 ■ 10 9 gauss called "atomic 
magnetic field", these sizes and energies are close to each other: ag = a#, E ~ uj L . We 
wish to study the spectrum of the hydrogen atom in magnetic fields much larger than B a . 
In this case the motion of the electron is mainly controlled by the magnetic field: it makes 
many oscillations in this field before it makes one in the Coulomb field of the nucleus. This 
is the condition for applicability of the adiabatic approximation, used for this problem for 
the first time in [T3] . 

The spectrum of a Dirac electron in a pure magnetic field is well known [15j; it admits a 
continuum of energy levels due to the free motion along the field: 

e 2 n = m 2 e + p 2 z + (2n + 1 + a z )eB , (22) 

where n = 0, 1, 2, a z = ±1 is the spin projection of the electron on z axis multiplied by 
two. For magnetic fields larger than B = m 2 /e, the electrons are relativistic with only one 
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exception: electrons belonging to the lowest Landau level (LLL, n = 0, a z = — 1) can be 
non-relativistic. 

In what follows we will study the spectrum of electrons from LLL in the Coulomb field 
of the proton modified by the superstrong B. The solution can be found in [TB] of the 
Schrodinger equation for an electron in a constant in time homogeneous magnetic field B in 
the gauge in which A = | B x r in cylindrical coordinates (p,z). The electron energies are: 



E 



p z n p ma z 



n + 



\m\ + m + 1 + a. 



eB pi 
h ~~~ 

m P 2m P 



(23) 



where n p = 0, 1, 2, ... is the number of nodal surfaces, m = 0, ±1, ±2, ... is the electron orbital 
momentum projection on the z axis (direction of the magnetic field) and a z — ±1. According 
to 1161. the LLL wave functions are: 



R 0m (p)= k24) 1+ H(| m |P 



p\ m \ e { im ^-p 2 l^ a H)) 



(24) 



P 



\p\, J \Ro m (p)\ 2 d 2 p = 1 , m = 0,-1,-2,.. 



We should now take into account the electric potential of the atomic nucleus located at 
p = z = 0. For an <C cib the adiabatic approximation can be used and the wave function 
should be looked for in the following form: 



^nOm(-l) = R()m(p)Xn(z) 



(25) 



where Xn{z) is the solution of a Schrodinger equation for an electron motion along the 
direction of the magnetic field: 



1 <f 



+ U eff {z) 



Xn{z) = E n Xn{z) 



2m e dz 2 ' ' ' 

Without screening the effective potential is given by the following formula: 



(26) 



Ueff(z) = -e 2 



\Rom(p)\' 

Vp^T^ 2 



d 2 P 



(27) 



which becomes the Coulomb potential for \z\ ^> a# 

Ueff{z) 



z 3> a H 



e 
\z\ 



(28) 



and is regular at z = 



Ueff(0) 



Off 



(29) 



To take screening into account we must use (fill to modify (1271) (see below). Since 
U e ff(z) = U e ff(—z), the wave functions are odd or even under reflection z — > —z; the 
ground states (for m = 0,-1, —2, ...) are described by even wave functions. 

In Fig. 5 the different scales important in the consideration of the hydrogen atom in 
strong magnetic field are shown. 
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analytical 
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Fig. 5. Landau radius an versus magnetic field B. 

To calculate the ground state of hydrogen atom in [T7j the shallow-well approximation 
is used: 



-2m e 



a B 



U{z)dz 



{m e ey2)ln 2 (B/(m 2 e e 3 )) 



(30) 



Let us derive this formula. The starting point is the one- dimensional Schrodinger equa- 
tion: 2 

- Y^ X{Z) + U &XW = E °x(z) (31) 
Neglecting Eq in comparison with U and integrating (|3TT) we get: 



X \a) = 2/i / U{z) X {z)dz 



(32) 



where we assume U(z) = U(—z), that is why x is even. 

The next assumptions are: 1. the finite range of the potential energy: U(z) ^ for 
a > z > —a; 2. x undergoes very small variations inside the well. Since outside the well 



x{z) ~ e V 2 v\ E °\ z } we readily obtain: 



\E \ = 2/i 



U{z)dz 



(33) 
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For 



/i\U\a 2 < 1 



(34) 



(condition for the potential to form a shallow well which means that the absolute value of 
the energy of ground level is much smaller than the absolute value of the potential in the 
well) we get that, indeed, \E Q \ <C \U\ and that the variation of x inside the well is small, 

Concerning the one-dimensional Coulomb potential, it satisfies this condition only for 
a <C l/(m e e 2 ) = an- This explains why the accuracy of log 2 formula ( 1301) is very poor. 

Much more accurate equation for atomic energies in strong magnetic field was derived by 
B.M.Karnakov and V.S.Popov [18] - It provides a several percent accuracy for the energies 
of EVEN states for H > 10 3 (H = B/(m 2 e e 3 )). 

Main idea is to integrate Shrodinger equation with effective potential from x = till 
x — z, where an « z « a B and to equate obtained expression for x'( z )/x{ z ) to the 
logarithmic derivative of Whittaker function - the solution of Shrodinger equation with 
Coulomb potential, which exponentially decreases at z » as- In this way in [18J the 
following equation was obtained: 

2 In (—) + In 2 - -0(1 + H) + 0(a H /z) = 
\a H J 

21n (^) + A + 21nA + 2V>(l--) + 4 7 + 2 In 2 + 0(z/a B ) , (35) 



E = -(m e e 4 /2)A 2 , 

where if>(x) is the logarithmic derivative of the gamma function. 
The energies of the ODD states are: 



E, 



odd 



A / 2 3 N 

m P e ^ / mze 
+ O 



2n 2 



B 



n 



1,2 



(36) 



(37) 



So, for superstrong magnetic fields B ~ m 2 /e 3 the deviations of odd states energies from 
the Balmer series are negligible. 

From (135]) we get an equation for A: 



ln(#) = A + 21nA + 2^(l 



A 



In 2 + 47 + 0(1 + M) 



(38) 



where ip(x) has simple poles at x = 0, —1, —2, .... So to reproduce large number at left hand 
side A should be large (ground level) or follow Balmer series (excited levels). 

When screening is taken into account an expression for effective potential transforms into 



U, 



2 I \Rom(p)\ 2 , 2 



6ml z _|_ e - A /(2/7r)e 3 B+6m2 z 



(39) 



Screening modifies the Coulomb potential at the distances \z\ < l/m e and since at these 
distances m e \U\a 2 = m e e 2 a < e 2 « 1, the approach leading to fl35l) still works. 

The modified Karnakov - Popov equation, which takes screening into account looks like: 



In 



H 



Sir / 



A + 21nA + 2^ (l - -A +ln2 + 47 + ^(1 + |m|) . 



(40) 
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We see that at H = B/{m 2 e e 3 ) ~ 3ir/e 6 ,B « 3iim 2 ./e 3 freezing of the energies occur: left 
hand side of (H0|) approach constant when B further grows. In particular, for a ground state 
at B » 3nml/e 3 we obtain: A = 11.2, E = -1.7 keV. 

Energy levels on which LLL is splitted in the hydrogen atom at B » 3irm 2 /e 3 are 
shown on Fig. 6. 
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-88- 
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Fig.6. Spectrum of hydrogen levels in the limit B » 3nm 2 /e z . Energies are given in 

rydberg units, Ry = 13.6 eV . 



5 Dirac equation with a screened Coulomb potential, 

Z = 1 

In the previous Section the spectrum of energies on which the lowest Landau level (LLL) splits 
in the proton electric field was found by solving the corresponding Schrodinger equation. 
Since the ground state energy of hydrogen in the limit of infinite B equals E Q = —1.7 keV, 
the use of the nonrelativistic Schrodinger equation is at least self consistent. However, the size 
an of the electron wave function for B > m^/e 3 in the direction transverse to the magnetic 
field is much smaller than the electron Compton wavelength, an = 1/y/eB < e/m e l/m e , 
which makes the nonrelativistic approach a bit suspicious (an = l/fn e for B = Bo). That is 
why in this Section we will study the ground state energy of the electron in a hydrogen-like ion 
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in the presence of an external magnetic field by analyzing the Dirac equation. Without taking 
screening into account this problem was considered in paper [19] (see also [20], were results 
obtained in [19] were reproduced), soon after it was found that a hydrogen-like ion becomes 
critical at Z « 170: the electron ground level sinks into the lower continuum (e < ~^ e ) 
and the vacuum becomes unstable by spontaneous e + e~ pairs production. These results 
were obtained by solving the Dirac equation for an electron moving in the field of a nucleus 
of finite radius. That the phenomenon of criticality can be studied only in the framework of 
the Dirac equation is an additional motivation for us to go from Schrodinger to Dirac. 

From the numerical solution of the Dirac equation for the ground electron level of a hy- 
drogen atom in the Coulomb potential we will find that the corrections to the nonrelativistic 
results are small and that the estimate SE = \E® — -Eq c1i | ~ (i?Q Ch ) 2 /m e works well. 

Let us parametrize bispinor which describes electron wave function in the following way: 

■ <P=(*) ■ X=( b A ■ (41) 



Substituting \1/ in the Dirac equation for the electron in an external electromagnetic field 
we obtain: 

+ (e + m — eip) [ ^ ] = 




"2 



Taking vector potential which describes constant magnetic field B directed along z axis 
in the form A = {—\By, \Bx, 0), we get: 





y + ix 
y — ix 

where p = \/x 2 + y 2 , 9 = arctan(y/x). Analogously we obtain: 

f> ( JL p-MJL - ie ' iB 9 

u 1 dz dp 

)d_ _|_ ie ie d 

dp p 89 dz 

Substituting two last expressions in the Dirac equation we get: 



(43) 




-i 



''"(-Mp + A + JA) ~m 




(45) 

Axial symmetry of electromagnetic field allows to determine 9 dependence of the functions 
c,- and hi. 




Cl (p,z) e^ M - l W e \ ( h\ _ f h(p,z) e^ 1 - 1 ' 2 ^ 
c 2 {p,z) e*(^+V2)e I. [b 2 ~[b 2 (p,z) e < M+1 ^ e 



(46) 



where M = ±1/2, ±3/2, ... is the projection of electron angular momentum on z axis. 
Substituting ( H61) in ( )45j) we get four linear equations for four unknown functions q and bi 
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(here and below C\ = ci(p, z), b\ = b\ (p, z) ...): 



/ , , M+ 1/2, e^p, , 

(e - m - eip)ci + i{-b lz - & 2p »2 = 

p 2 

/ , M-l/2, eEp, , . 

(e - m - e<p)c 2 + H-Oip H h + —r-h + b 2z ) = 

p 2 

M+l/2 eBp . n 

(e + m- etp)bi + i{-c lz - c 2p c 2 — c 2 ) = 

p z 

x, M- 1/2 e£p 

(e + m - e<p)6 2 + i(-c lp + '—c x + — ^c x + c 22 ) = , (47) 

p z 

where 6 lz = db%/dz, b\ p = dbi/dp, ... Ground energy state has s z = — 1/2, l z = 0. Taking 
M = — 1/2 we should look for solution of (|47|) with c x = 61 = 0: 

' b 2p + e -fb 2 = 

(48) 

, c 2p +^c 2 = , 

f (e - to - ey?)c 2 + 162* = , 4Q , 
[ (e + m- etp)b 2 + ic 2z = . 

The dependence on p is determined by (j4"8~l) : 

f 6 2 (p,z) = e-^ 2 /4 ( _, )/W 

1 c 2 (p,z) = e-^ (z) . (50) 



Substituting the last expressions in (T4l?j) and averaging over fast motion in transverse 
to the magnetic field plane we obtain two first order differential equations which describes 
electron motion along magnetic field in an effective potential V(z): 

g z - (e + m e - V)f = 0, , , 

f z + (e- m e - V)g = 0, l ° > 



Ze 2 ? ex P 



2ai 



9{2) = -4; J ^FrW pdp ■ (52) 

At large distances \z\ ^> a# the effective potential equals Coulomb, and the solutions 
of the equations (15T!) exponentially decreasing at \z\ — > 00 are linear combinations of the 
Whittaker functions. At short distances the equations (I5TI) can be easily integrated for 
|t^(z)| 3> |e ± m e \ (as far as |e| < m e condition for IV'(z)! will be for sure valid for |V(2)| > 
2m e , which is equivalent to the following inequality: z <C zTe 2 /(2m e )), where they looks like: 

9z + Vf = , f z - Vg = . (53) 

The result of the integration is: 

g(z) = Bi cos w(z) + B 2 sin w(z) , 

f(z) = B\ sin w(z) — B 2 cosw(z) (54) 
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where 

w( z ) = I V(z)dz (55) 



o 

and Bi, B2 are normalization constants. 

The functions g(z) and f(z) have opposite parities; for the ground state g(z) should be 
even, so B% = 0, and matching logarithmic derivatives at the point zq we obtain: 

- V(z ) tan w(z ) = -jj- In g(z ) , (56) 

Ze 2 /(2m e ) > z > a H (57) 

Substituting proper combination of the Whittaker functions for g(z) we obtain an algebraic 
equation for the ground state energy (it coincides with Eq. (22) in [19j in the limit R/clh *C 1, 
where R is the nucleus radius): 



Zi'\ u \ 2 v ' m !_ £ - | +argT ( ; = + iZe 2 \ 

feB } \ Jm*-e 2 



(fry, _|_ £ \ 7T 
J— - argT(l + 2iZe 2 ) (In 2 + 7) = - + mr , (58) 
V m e — e J 2 2 

where 7 = 0.5772... is the Euler constant, and the argument of the gamma function is given 
by 

arg T(x + iy) = -jy + X] (f _ arctan - + ^ - - j . (59) 

For the ground level at e > one should take n = 0, while for e < it should be changed to 
n = —1. 

According to (!58j) when the magnetic field increases the ground state energy goes down 
and reaches the lower continuum. 

A matching point exists only if B ^> Am 2 / (e(Ze 2 ) 2 ) (see ([57)1 ) and (1581) is valid only for 
these values of the magnetic field. 

Thus, without taking screening into account, from fl58|) we can obtain the dependence 
of the ground state energy of a hydrogen atom on the magnetic field for B ^> Am 2 /e 5 . 
Screening modifies the Coulomb potential at distances smaller than the electron Compton 
wave length, and from the condition |V(l/m e )| ^> 2m e we get Ze 2 » 2. It means that at 
B > 3nm 2 /e 3 the phenomena of screening does not allow to find analytically the ground state 
energy. In order to find the ground state energy at B < Am 2 / e 5 and to take screening into 
account the equations (!5T!) were solved numerically. This system can be transformed into 

1 /2 

one second order differential equation for g(z). By substituting g(z) = (e + m e — V^j x( z ) 
a Schrodinger-like equation for the function x( z ) was obtained in [T§]rl 

^ + 2m e (E-U) X = , (60) 



lr This trick was exploited by V.S. Popov for the qualitative analysis of the phenomenon of critical charge 
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E = fl^I , „ = ±- 9 - —V 2 + ? - + 3/8(f " )2 



2m e ' m e 2m e Am e (e + m e — V) m e (e + m e — V^) 2 ' 

where £ is the energy eigenvalue of the Dirac equation and V(z) is given in ( l52l) . An equation 
(|60l) was integrated numerically. Let us note that, while for z ^> l/m e the last three terms 
in the expression for U are much smaller than the first one (the only one remaining in the 
nonrelativistic approximation), at z < l/m e the relativistic terms dominate and are very big 
for B ^> Bq at z ~ a# which makes numerical calculations very complicated. 

In Tabled] the results for the ground state energy of a hydrogen atom without screening 
are presented. The values of the magnetic field in units of -B^are given in the first column 
while in columns 2-5 the values of A are given. By definition 



o 
2 



E = £ -^- = -^A 2 . (61) 
2m e 2 v ' 



From Tabled] we see that: 



1. the results of the numerical integrations of the Schrodinger and Dirac equations coincide 
within four digits: 

• with the analytical Karnakov-Popov formula for the ground state energy (n p = 
m = 0) [18J in the case of the Schrodinger equation; 

• with formula fl58|) for Z = 1 in the case of the Dirac equation; 

2. for the relativistic shift of energy the following estimate works: 

-^Dirac — -Kschr ~ -^Schr— , (62) 

5A~e 4 A 3 /4 . 

To take screening into account, the following formula for the effective potential should 
be used in (160]) instead of ( 152]) : 



2 r , , n °? „-p 2 /2a 2 . 



a 



H 



1 _ e-v^fN + e -VVM* 3B +e™ 2 M f e " v " odp , (63) 

- J VP + z 



where Z = 1 for hydrogen. 

The freezing of the ground state energy is due to a weaker singularity of the potential 
with screening (163 j) at z — > for B — > oo than that of the potential without screening 
(I52p . While the non-screened potential behaves like 1/z at small z, the screened potential is 
proportional to 5(z) because, when B — > oo, the width of the region where it behaves like 
1/z shrinks to zero [2] . 

In Table [2] the results of the analytical formula for A with the account of screening for 
the Schrodinger equation ( HO]) are compared with the results of the numerical integration of 
the Dirac equation. We see that in the case of screening the relativistic shift of energy is 
also very small, and due to it the ground state energies become a little bit higher, just like 
without taking screening into account. The freezing of the ground state energy occurs at 
B/B = 10 3 + 10 4 , when B w 3irm 2 Je 3 . 



2 Let us note that the definition of A used in [19] differs from our: aQ2] 
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Table 1: Values of A for Z = 1 without screening obtained from the Schrodinger and Dirac 
equations. They start to differ substantially at enormous values of the magnetic field. 



B/B 


KP-equation 


Numerical results 


T?n /liro[\ 

bq. (|58j) 


Numerical results 




^ocnroaingerj 


(Schrodinger) 


(Dirac) 


(Dirac) 


i nO 
ill 


5.737 


C TO C 

5.735 


5.735 


5.734 


10 


7.374 


7.374 


T o Tn 

7.370 


T O T1 

7.371 


10 


A 1 /II 

9.141 


y.141 


n 1 op 

9.136 


n 1 o c 

9.135 


ill 


11.00 


1 1 nn 

11.00 


1 n nn 

10.99 


i n nn 

10.99 


10 


12.93 


1 o no 

12.93 


1 o n 1 

12.91 


1 o n 1 

12.91 


10 5 


14.91 


14.91 


14.88 


14.88 


10 6 


16.93 


16.93 


16.89 


16.89 


10 7 


18.98 


18.98 


18.93 


18.92 


10 8 


21.06 


21.05 


20.98 


20.98 


10 9 


23.16 


23.15 


23.05 


23.05 


10 10 


25.27 


25.27 


25.14 


25.13 


10 11 


27.40 


27.40 


27.23 




10 12 


29.54 


29.54 


29.33 




10 15 


36.03 




35.64 




10 20 


46.99 




46.11 




10 25 


58.07 




56.40 




10 30 


69.22 




66.38 




10 35 


80.43 




75.98 




10 40 


91.67 




85.10 




10 45 


102.95 




93.67 




1Q 50 


114.25 




101.62 




10 55 


125.57 




108.89 





6 Screening versus critical nucleus charge 

According to [TH] nuclei with Z > 40 become critical in an external B (for smaller Z the 
values of at which the criticality is reached become smaller than the nucleus radius, the 
Coulomb potential diminishes and thus the ground level does not reach the lower continuum). 

In Table [3] one can see the dependence of the ground state electron energy e on the ex- 
ternal magnetic field for Z = 40. The numerical solutions of (167J1) are in good correspondence 
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Table 2: Values of A for Z — 1 with screening. 



B/B 


hiq. (1401) 


Numerical results 


Numerical results 




^ocnroctmgerj 


^ocnroamgerj 


(Dirac) 


i nO 
10 


5.7 


5.7 


5.7 


1 nl 
10 


7.4 


7.4 


7.4 


in 2 


Q 1 


Q 1 


Q 1 
y. i 


10 3 


10.5 


10.6 


10.6 


10 4 


11.1 


11.2 


11.2 


10 5 


11.2 


11.3 


11.3 


10 6 


11.2 


11.4 


11.3 


10 7 


11.2 


11.4 


11.3 


10 8 


11.2 


11.4 


11.3 



Table 3: Values of So/m e for Z = 40. 



B/B 


Eq. §8j) 


Numerical results 


Numerical results 




(Dirac) 


(Dirac) 


with screening (Dirac) 


10° 


0.819 


0.850 


0.850 


10 1 


0.653 


0.667 


0.667 


10 2 


0.336 


0.339 


0.346 


10 3 


-0.158 


-0.159 


-0.0765 


10 4 


-0.758 


-0.759 


-0.376 


2 • 10 4 


-0.926 


-0.927 


-0.423 




at B/B f 


a 2.85 ■ 10 4 , e = -m e 




10 5 






-0.488 


10 6 






-0.524 


10 7 






-0.535 


10 8 






-0.538 



Table 4: Values of freezing ground state energies for different Z from the Schrodinger and 
the Dirac equations. In order to find the freezing energies we must take B/Bq 3> 37r/e 2 . In 
numerical calculations we took B/B = 10 8 . 



z 


/ numerical 

(E f A , keV 

V U JSr.hr 


/ f r \ numerical 

( En — m e ) , keV 

V / Dirac 


1 


-1.7 


-1.7 


10 


-88 


-87 


20 


-288 


-273 


30 


-582 


-519 


40 


-966 


-787 


49 




-1003 



with the values of Eq obtained from fl58|) . The numerical results with screening are shown in 
the last column; we see that freezing occurs in the relativistic domain e ps — m e /2 and the 
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ground level never reaches lower continuum, Eq > —m e . 

In Table H] we compare freezing energies for different Z obtained numerically from the 
nonrelativistic Schrodinger equation and from the Dirac equation. We see that for Z > 20 
the freezing occurs in the relativistic regime, where the Schrodinger equation should not be 
used. Let us stress that the value of B at which the freezing occurs does not depend on Z. 

From (15 8p we obtain in the limiting case e — > — m e an equation which defines the value 
of the magnetic field at which a nucleus with charge Z becomes critical without taking 
screening into account (it coincides with Eq. (32) from [T9]): 

B ofy 2 , 2 / , 7r-2argr(l + 2^ cr e 2 ) ^ 

- = 2(Z cr e ) exp |^- 7 + — j . (64) 

This equation is used to calculate the numbers in the second column of Table EJ 



Table 5: Values of B/Bq at which Eq = —m e according to the Dirac equation and nuclei be- 
come supercritical without (column 2,3) and with (column 4) taking screening into account. 



7 


Eq. ([MD 


Numerical results 


Numerical results 






without screening 


with screening 


90 


118 


116 


122 


85 


157 


154 


164 


80 


213 


210 


229 


75 


301 


297 


335 


70 


444 


438 


527 


65 


689 


681 


923 


60 


1144 


1133 


1964 


55 


2068 


2053 


6830 


54 


2357 


2340 


10172 


53 


2699 


2681 


17012 


52 


3107 


3087 


35135 


51 


3594 


3572 


1.20 ■ 10 5 


50 


4181 


4157 


1.14 ■ 10 7 


45 


9826 


9787 




40 


28478 


28408 




35 


1.12 • 10 5 


1.12 • 10 5 




30 


6.99 • 10 5 


6.98 • 10 5 




25 


9.27- 10 6 


9.27 • 10 6 





From Table [5] we see that with the account of screening only the atoms with Z > 52 
become supercritical at the values of B/Bq shown in the fourth column. Because of screening 
a larger B is needed for a nucleus to become supercritical and the nuclei with Z < 52 never 
reach supercriticality. This phenomenon is illustrated in Fig. 7. 
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10 a 
10 7 
10 6 

B/B u io 5 

10 4 



10 3 



iu 20 30 40 50 60 70 80 90 

Fig. 7. The values of Bj? r : a) without screening according to (64\ ), dashed (green) line; b) 
numerical results with screening, solid (blue) line. The dotted (black) line corresponds to 
the field at which an becomes smaller than the size of the nucleus. 

From Tables [I] El and [5] we see that (1581) is very good in describing the dependence of 
the energy on the magnetic field; at least a numerical integration produces almost identical 
results. In Table [6] we demonstrate several cases where the accuracy of (I58p is not that 
good. It happens at low B/Bq since the matching condition B > 4m 2 /(e(Ze 2 ) 2 ) fails and 
when Eq is relativistic. However, B should not be too low to make the adiabaticity condition 
clb 3> an, or B 3> (Ze 2 ) 2 m 2 e / e applicable. 



Table 6: Values of e /m e at B/B = 5. 



z 


Eq. (EHD 


Numerical results 




(Dirac) 


(Dirac) 


90 


0.2050 


0.2512 


80 


0.3096 


0.3539 


70 


0.4139 


0.4542 


60 


0.5171 


0.5516 


50 


0.6185 


0.6454 


40 


0.7165 


0.7349 


30 


0.8086 


0.8188 


20 


0.8914 


0.8952 


10 


0.9596 


0.9601 


1 


0.998745 


0.998745 



Textbooks [22] contain detailed consideration of the phenomenon of critical charge. 
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7 Conclusions 



An analytical formula for the Coulomb potential 3>(z) in a superstrong magnetic field has 
been derived. It reproduces the results of the numerical calculations made in [2] with good 
accuracy. Using it, an algebraic formula for the energy spectrum of the levels of a hydrogen 
atom originating from the lowest Landau level in a superstrong B has been obtained. The 
energies start to deviate from those obtained without taking the screening of the Coulomb 
potential into account at B > 3nml/e 3 ~ 6 ■ 10 16 gauss and the energy of ground state in 
the limit B — > oo remains finite. 

A magnetic field plays a double role in the critical charge phenomenon. By squeezing 
the electron wave function and putting it in the domain of a stronger Coulomb potential it 
diminishes the value of the critical charge substantially [19] . However, for nuclei with Z < 52 
to become critical such a strong B is needed that the screening of the Coulomb potential 
occurs and acts in the opposite direction: the electron ground state energy freezes and the 
nucleus remains subcritical in spite of growing B. 

I am very grateful to Bruno Machet and Sergey Godunov for fruitful collaboration. 
I am partially supported by the RFBR grants 11-02-00441 and 12-02-00193, by the grant 
of the Russian Federation government 11.G34.31.0047, and by the grant NSh-3172.2012.2. 
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